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Abstract 


The electrical performance characteristics of 2 and 10 fl-cm N/P-type silicon 
solar cells were measured at simulated solar intensities of 5, 50, 100, 140, 250, 400, 
550, 700, and 850 mW/cm 2 . At each intensity, the temperature was varied in incre- 
ments of 20° between extremes of + 160 and — 160°C. Short-circuit current, open- 
circuit voltage, and maximum power are presented in graphical format. Also 
described are three methods for predicting solar cell electrical performance as a 
function of temperature and intensity. Two of the methods are suitable for use 
at extreme temperature-intensity conditions. These methods were used success- 
fully to predict the performance of the I 8C — V oc transducer on board the Mariner 
Mars 1971 spacecraft. 
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Parametric Performance Characteristics and Treatment 
of Temperature Coefficients of Silicon 
Solar Cells for Space Application 


I. Introduction 

Silicon solar cell arrays have proven themselves as the 
most practical, reliable, and efficient source of electrical 
power for planetary spacecraft missions flown to date. One 
of the important prerequisites for successful design and 
operation of a solar array for a particular mission is the 
ability to accurately predict its performance character- 
istics over the appropriate solar intensity and temperature 
ranges. A major problem that limits prediction accuracy, 
particularly for missions involving temperature-intensity 
extremes, is the lack of good experimental test programs 
to provide statistically significant empirical test data. This 
lack of test data has limited design tradeoff studies and 
the development of solar array technology. In addition, 
methods for predicting solar cell performance for near- 
Earth missions are not suitable for application to extreme 


temperature-intensity mission predictions without sacri- 
ficing accuracy. As a result, JPL has, over the years, en- 
gaged in a comprehensive study of the electrical and 
mechanical performance behavior of both conventional 
and nonconventional silicon solar cells for missions over 
the appropriate range of temperature-intensity conditions 
from Mercury to Jupiter. Also, JPL has developed several 
suitable methods for predicting solar cell performance 
under extreme temperature-intensity conditions. The 
methods, which are well suited for computer computa- 
tion, were developed with a minimum use of derived 
quantities; the identity of the raw parametric data on 
which they are based was therefore preserved. A con- 
densed version of this report (Ref. 1) was presented at the 
7th IECEC held in San Diego, California, on Septem- 
ber 25-29, 1972. 
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II. Experimental Approach 

A. Solar Cell Description 

The silicon solar cells of 2 and lOfl-cm base resistivity 
discussed in this report were fabricated by Centralab and 
Heliotek, respectively. In order to maintain a common 
base bne, the cells are similar in all other respects. Both 
cell types are N/P, are typically 0.046 cm (18 mils) thick, 
have a surface area of 4 cm 2 , have solder-coated titanium- 
silver vacuum-evaporated bar contacts, and were assem- 
bled using a solder reflow technique (Ref. 2). A fused silica 
coverglass, Comings glass works 7940, 0.051 cm (20 mils) 
thick, was bonded to each cell with General Electric RTV- 
602/SRC05 silicone adhesive. 

B. Apparatus 

The electrical parametric studies performed consisted 
of generating current voltage characteristics of two basic 
types of silicon solar cells over the appropriate range of 


temperature-intensity, equivalent to that encountered 
during Mercury and Jupiter missions. Figure 1 shows an 
overall view of the test facility used throughout the course 
of these investigations. A multi-cell thermal vacuum test 
chamber was employed in conjunction with a solar simu- 
lator, which approximates the intensity and spectrum of 
space sunlight. 

A Spectrolab Model X25 Mark II, close-filtered, dry 
nitrogen-purged solar simulator was used throughout the 
test program. In order to achieve the broad range of solar 
intensities required, two separate lenticular optical units 
were used, one for solar intensities from 250 to 5 mW/cm 2 
and the other for intensities from 250 to 850 mW/cm 2 . 
Intermediate intensities were obtained by interposing vari- 
ous types of screens within the light path and adjusting the 
current to the xenon source. Balloon flight standard cells 
comprising an intensity reference cell and two separate 
cells having narrow-bandpass filters were mounted within 
the test chamber. The latter cells were used to determine 



Fig. 1. Solar energy developmental and test facility 
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the spectral quality of the solar simulator by observing its 
red-to-blue ratio during test. 

Electrical connection to each cell was made with four 
wires in order to prevent voltage drop across the current- 
carrying pair from influencing the voltage readings. A 
typical solar cell parametric test plate composed of fil- 
tered cells, thermocouples, interconnectors, printed circuit 
board, and wiring harness is shown in Fig. 2 as viewed 
through the quartz window of the test chamber. The solar 
cells are bonded to a 0.318-cm (0.125-in.) thick copper 
plate using General Electric RTV-560 silicone adhesive 
and primer. The electrical characteristics of each cell 
tested are obtained in the form of a current-voltage curve 
using a Hewlett Packard Model 7030 XY recorder. The 
short-circuit current and open-circuit voltage parameters 
are obtained on a five-place readout, integrating-type 
digital voltmeter. The operating temperature of the 
solar cells was controlled within ±0.5°C using copper- 
constantan thermocouples and a proportional temperature 
controller. Because of the extreme temperature range 
covered in this study ( — 160 to 4-160°C), solar cells were 
soldered to tinned Kovar interconnectors. (Kovar is a metal 
alloy which approaches the thermal expansion coefficient 
of silicon.) The cell temperature was varied in increments 
of 20° between temperature extremes of 4-160 and 
— 160°C. 

III. Graphical Presentation of Parametric Data 

All data presented in this report are based on the aver- 
age of 15 cells in each of two groups, one from Heliotek 
and the second from Centralab. Figure 3 shows the results 
of plotting short-circuit current (7 SC ) vs. solar intensity, 
with temperature as a parameter, for both 2 and 10 fl-cm 
cells. The 2 O-cm cell l sc is linear with intensity over the 
temperature-intensity range studied. However, the 10 
O-cm cell I sc does not behave linearly above 250 mW /cm 2 
and 100°C. The nonlinear behavior of the 10 O-cm cells is 
probably due to the higher series resistance of these cells. 
The linear behavior of the 2 O-cm cell makes it a much 
better candidate for a standard cell used in adjusting and 
monitoring solar simulator intensity. Also shown in Fig. 3 
is the I ac dependence on temperature. At least three 
separate factors contribute to reduced l sc at lower tem- 
perature. As temperature decreases, the band gap increases 
and fewer incident photons can generate current. Dash 
and Newman (Ref. 3) have shown that light absorption 
(necessary for the generation of current) decreases with 
temperature from 27 to — 196°C. Also, reduced collection 
of current carriers as temperature decreases results from 
a decrease in diffusion lengths at lower temperature. 



Fig. 2. Test specimens and standard cells as viewed through 
the quartz window of the test chamber 


Figure 4 shows the results of plotting open-circuit 
voltage V oc vs. the log 10 of solar intensity, with tempera- 
ture as a parameter, for both 2 and 10 fl-cm cells. The V oc 
of both cell types is a logarithmic function of intensity 
over the range of 25 to 850 mW/cm 2 . 

Maximum power P mp of both 2 and 10 fl-cm cell types 
is plotted vs. solar intensity, with temperature as a param- 
eter, in Fig. 5. As expected, P mp increases with intensity, 
since both I sc and V oc increase with intensity. However, 
at higher intensities (> 250 mW/cm 2 ), P mp of the 10 fl-cm 
cells falls off with respect to P mp of the 2 fl-cm cells. Thus 
the 10 fl-cm cells would not be good candidates for near- 
Sun missions. Also seen in Fig. 5 is the increase in P mp 
with decreased temperature. Thus, the rate of increase in 
V oc with decreased temperature is dominant over the rate 
of decrease of other parameters such as The spread in 
maximum power, as reflected by standard deviations, is 
seen to be greater at the low and high intensity extremes. 
This is as expected, since the cells used in this study were 
selected such that their maximum powers were matched 
at 1 AU equivalent solar intensity and 28° C temperature. 
The variability exhibited in P mp would cause serious mis- 
match problems in an array, which would result in a 
decrease in operating efficiency. As planetary exploration 
missions are extended toward inner and outer planets, 
further attention must be focused on cell fabrication 
processes, proper selection of materials, and cell screening 
techniques in order to ensure optimum performance. 
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IV. Solar Cell Modeling for Performance 
Predictions 

The parametric presentation of the solar cell electrical 
performance data in Figs. 3 through 5 is useful in terms of 
describing both a quantitative and qualitative relationship 
between solar cell output and temperature-intensity con- 
ditions. However, the quantitative description is limited in 
resolution in accordance with the number of temperature- 
intensity test conditions examined. Thus, in order to pre- 
dict the performance characteristics of silicon solar cells 
over a continuous range of temperature-intensity condi- 
tions, it is necessary to model the electrical performance 
of the solar cells. Three methods for modeling the para- 
metric data in order to predict solar cell performance are 
discussed below, together with the characteristic advan- 
tages and disadvantages of each. In general, the selection 
of one method over another involves trading off accuracy 
with complexity. 

A. Method 1: Analytical Expression 1 

This method has been previously described by Sand- 
strom (Ref. 4). The following equations are used to 
translate short-circuit current and open-circuit voltage 
from a reference temperature-intensity T 0 ,H 0 to a new 
temperature-intensity T it Hf. 


ho (T u Hi) = I, c (H 0 , T 0 ) + «„< (Ti - T 0 ) (1) 

V oc (Ti, Hi) = V oc (H 0 , To) - (Ti - To) - M SC R, 

(2) 

A/ sc = I sc (Ti, Hi) - Isa (To, H 0 ) (3) 

The series resistance R g is treated as a constant. The 
short-circuit current temperature coefficient a and the 
open-circuit voltage temperature coefficient /? are defined 
as the magnitude of the slope of the curves which express 
I sc and V oc , respectively, as a function of temperature at a 
given intensity. That is. 


di sc 

( 4 ) 

aHi = ~dJ 

dV 0 c 

( 5 ) 

Pm ~ dT 
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Fig. 5. Maximum power output as a function of 
intensity and cell temperature 


at constant intensity Hi. 

Although this method is easy to apply, it suffers several 
shortcomings. The coefficients a and /? do vary with tem- 


perature, and thus the valid application of Eqs. (1) and (2) 
is limited to a relatively small temperature range. Also, 
Eq. (2) does not explicitly reveal the logarithmic depen- 
dence of V oc on intensity which was shown in Fig. 4. 
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Table 1. Average short-circuit currents and standard deviations for 2 fi-cm cells 


Tempera- 





Solar intensity,' 

1 mW/cm 2 





hire, °C 

5.00 

25.00 

50.00 

100.00 

140.00 

250.00 

400.00 

550.00 

700.00 

850.00 

-160 

3.62 











(0.28) 

- 

- 

- 

- 

- 

- 

- 

- 

- 

-140 

3.67 

17.72 










(0.27) 

(1.05) 

- 

- 

- 

- 

- 

- 

- 

- 

-120 

3.53 

17.87 

36.69 









(0.28) 

(0.95) 

(2.06) 

- 

- 

- 

- 

- 

- 

- 

-100 

3.62 

18.29 

37.63 

73.79 








(0.27) 

(0.89) 

(1.90) 

(3.08) 

- 

- 

- 

- 

- 

- 

-80 

3.77 

19.18 

38.75 

77.18 

109.62 







(0.27) 

(0.73) 

(1.48) 

(2.66) 

(4.41) 

- 

- 

- 

- 

- 

-60 

3.95. 

19.81 

40.48 

81.93 

116.12 

209.25 






(0.26) 

(0.60) 

(1.20) 

(1.96) 

(3.07) 

(6.41) 

- 

- 

- 

- 

-40 

4.13 

20.79 

42.16 

85.51 

121.98 

221.10 






(0.26) 

(0.46) 

(1.03) 

(1.53) 

(2.68) 

(5.46) 

- 

- 

- 

- 

-20 

4.33 

21.25 

43.16 

88.30 

126.05 

228.45 






(0.26) 

(0.73) 

(0.87) 

(1.74) 

(2.33) 

(5.24) 

- 

- 

- 

- 

0 

4.41 

22.39 

43.91 

91.96 

129.33 

234.35 






(0.22) 

(0.35) 

(0.77) 

(1.33) 

(2.28) 

(5.38) 

- 

- 

- 

- 

20 

4.46 

23.23 

45.39 

93.45 

133.09 

239.20 






(0.27) 

(0.34) 

(0.93) 

(1.02) 

(2.10) 

(5.14) 

- 

- 

- 

- 

40 

4.86 

23.55 

46.51 

94.83 

134.85 

239.11 






(0.17) 

(0.34) 

(0.84) 

(1.15) 

(1.96) 

(4.81) 

- 

- 

- 

- 

60 

4.90 

23.83 

47.17 

95.77 

134.16 

242.15 

374.42 

514.25 

672.31 

835.33 


(0.16) 

(0.34) 

(0.74) 

(1.27) 

(1.82) 

(4.70) 

(9.93) 

(14.81) 

(11-14) 

(13.24) 

80 





136.29 

245.89 

377.79 

521.71 

672.77 

842.89 


— 

- 

— 

— 

(1.96) 

(4.53) 

(9.15) 

(11.67) 

(10.16) 

(12.22) 

100 





136.02 

249.07 

378.44 

525.64 

676.17 

850.05 


- 

- 

- 

- 

(1.98) 

(4.57) 

(10.51) 

(11.30) 

(9.50) 

(12.27) 

120 





137.28 

249.62 

387.15 

526.78 

683.42 

856.10 


— 

- ■ 

- 

- 

(1.96) 

(4.43) 

(10.17) 

(8.41) 

(9.70) 

(11.94) 

140 





139.13 

252.56 

395.12 

537.79 

679.29 

857.11 


— 

— 

- 

- 

(2.03) 

(4.59) 

(10.14) 

(11.80) 

(10.37) 

(13.94) 

160 






255.21 

398.19 

532.66 

685.62 

860.11 


- 

- 

- 

- 

- 

(4.29) 

(10.06) 

(8.22) 

(12.54) 

(19.39) 

Standard deviations are given in parentheses. 


Equations (1) and (2) have been used to successfully 
predict solar cell performance for near- Earth conditions; 
however, because of the limitations discussed above, 
methods 2 and 3 were developed. 

B. Method 2: Analytical Expression 2 

This method was developed for the purpose of provid- 
ing an analytical expression well suited for computer com- 
putation which could be used for predicting I sc and V oc 
over a wide range of both temperature and intensity. Also, 
it was desired to take into account the relative uncertainty 
of the parametric data. The I sc and V oc parametric data 
(average output of 15 cells with standard deviations) of 
the 2 and 10 Q-cm cells are tabulated in Tables 1 through 4. 


The V oc and the I ac models were developed in a similar 
manner. The development of the V oc model will be dis- 
cussed in detail. To avoid repetition, only the final results 
of the I sc model will be described. 

1. Development of V oc Model. The first step was to 
obtain an expression which described V oc performance as 
a function of intensity at constant temperature T { . The 
following equation was found to fit both the 2 and 10 Q-cm 
parametric data quite well over the intensity range of 25 
to 850 mW/cm 2 : 

V oc (T i; H) = A Ti + B t . logio H (6) 
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Table 2. Average short-circuit currents and standard deviations for 10 £l-cm cells 


Tempera- 




Solar intensity, a 

mW/cm 2 





ture, °C 

5.00 

25,00 

50.00 

100.00 

140.00 

250.00 

400.00 

550.00 

700.00 

850.00 

-160 

4.38 











(0.22) 

- 

- 

- 

- 

- 

- 

- 

- 

— 

-140 

4.48 

21.49 










(0.19) 

(0.63) 

- 

- 

- 

- 

- 

- 

- 

- 

-120 

4.62 

22.04 

44.30 









(0.16) 

(0.44) 

(0.63) 

- 

- 

- 

- 

- 

- 

- 

-100 

4.73 

22,52 

45.30 

91.05 








(0.14) 

(0.34) 

(0.46) 

(1.42) 

- 

- 

- 

- 

- 

- 

-80 

4.79 

22.88 

45.99 

93.01 

130.65 







(0.13) 

(0.34) 

(0.36) 

(1.35) 

(2.49) 

- 

- 

- 

- 

- 

-60 

4.87 

23.52 

46.29 

95.81 

132.65 

237.40 






(0.14) 

(0.34) 

(0.31) 

(1.00) 

(2.50) 

(1.78) 

- 

- 

- 

- 

-40 

4.73 

23.90 

47.33 

96.92 

135.39 

240.88 






(0.11) 

(0.34) 

(0.30) 

(0.95) 

(3.29) 

(1.95) 

- 

- 

- 

- 

-20 

4.73 

24.36 

48.50 

98.83 

137.41 

245.55 






(0.12) 

(0.36) 

(0.27) 

(0.99) 

(2.96) 

(2.14) 

- 

- 

- 

- 

0 

4.81 

24.79 

48.79 

100.55 

141.05 

252.99 






(0.13) 

(0.32) 

(0.33) 

(1-21) 

(2.91) 

(2.15) 

- 

- 

- 

- 

20 

4.79 

25.26 

50.43 

102.64 

145.57 

258.54 






(0.13) 

(0.34) 

(0.38) 

(1.13) 

(3.16) 

(1.82) 

- 

- 

- 

- 

40 

4.89 

25.64 

51.46 

104.05 

145.39 

261.90 






(0.14) 

(0.36) 

(0.41) 

(1.12) 

(3.20) 

(2.08) 

- 

- 

- 

- 

60 

4.98 

25.92 

52.25 

105.31 

147.02 

265.02 

412.71 

575.35 

717.70 

850.13 


(0.13) 

(0.36) 

(0.39) 

(1.15) 

(3.26) 

(2.13) 

(4.03) 

(6.31) 

(9.17) 

(13.95) 

80 





148.34 

267.92 

415.52 

577.60 

717.05 

840.63 


- 

- 

- 

- 

(3.57) 

(2.22) 

(4.17) 

(8.04) 

(9.63) 

(19.48) 

100 





149.85 

265.04 

416.19 

576.09 

708.15 

820.55 


- 

- 

- 

- 

(3.50) 

(1.96) 

(4.48) 

(8.20) 

(11.87) 

(27.22) 

120 





154.22 

266.89 

415.55 

567.68 

686.22 

779.98 


- 

- 

- 

- 

(3.53) 

(3.30) 

(5.04) 

(11.13) 

(17.59) 

(33.46) 

140 





154.95 

274.97 

411.10 

538.27 

642.26 

723.32 


- 

- 

- 

- 

(3.46) 

(2.15) 

(7.02) 

(15.03) 

(24.80) 

(39.31) 

160 






257.56 

387.45 

496.61 

580.27 

646.08 


- 

- 

- 

- 

- 

(4.98) 

(11.49) 

(21.22) 

(30.41) 

(39.31) 

“Standard deviations are given in parentheses. 


The coefficients A T . and B r . were found for each tem- 
perature by a least-squares fit to the parametric data, 
weighted according to the inverse standard deviation u of 
each data point. That is, at temperature 7%, A T . and B T< 
were found such as to minimize the following summation: 

f v. e (T < ,H,)-(A r< + B r< logH,) 

( <r(Ti,Hj) 

No attempt was made to counterweight the least-squares 
fit to “neutralize” the weighting incurred by taking the log 
of intensity. 




The temperature dependence of the A T . and the B r . was 
found empirically. Polynomials of first through fifth order 
were fitted individually to both the A r . and B T . vs. tem- 
perature plots. Fifth-order polynomials gave the best fit 
for both A(T) and B(T). Thus, A(T) and B(T) take the 
following form: 


A (T) = o 0 + aj' + aj* + a 3 T 3 + aj * + a 5 T 5 (8) 

B ( T ) = b 0 + b,T + b 2 T* + b 3 T 3 + bj* + h 5 T 5 (9) 
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Table 3. Average open-circuit voltages and standard deviations for 2 S!-cm cells 


Tempera- 





Solar intensity, 1 

1 mW/cm 2 





ture, °C 

5.00 

25.00 

50.00 

100.00 

140.00 

250.00 

400.00 

550.00 

700.00 

850.00 

-160 

834.48 











(46.90) 

- 

- 

- 

- 

- 

- 

- 

- 

- 

-140 

801.05 

866.62 










(39.67) 

(23.51) 

- 

- 

- 

- 

- 

- 

- 

- 

-120 

772.26 

837.71 

849.79 









(32.50) 

(13.47) 

(15.58) 

- 

- 

- 

- 

- 

- 

- 

-100 

735.09 

805.89 

824.37 

830.75 








(23.07) 

(7.01) 

(8.40) 

(10.42) 

- 

- 

- 

- 

- 

- 

-80 

693.93 

765.08 

785.50 

800.18 

808.45 







(18.11) 

(5.00) 

(5.30) 

(5.05) 

(5.99) 

- 

- 

- 

- 

- 

-60 

653.18 

725.24 

746.55 

765.28 

771.85 

786.66 






(14.75) 

(3.98) 

(3.31) 

(4.07) 

(4.51) 

(6.75) 

- 

- 

- 

- 

-40 

608.94 

681.66 

705.83 

724.30 

737.19 

752.17 






(13.00) 

(3.11) 

(2.53) 

(3.00) 

(3.68) 

(6.10) 

- 

- 

- 

- 

-20 

565.54 

537.40 

661.98 

682.20 

696.16 

713.61 






(10.81) 

(2.65) 

(2.11) 

(2.28) 

(2.72) 

(4.32) 

- 

- 

- 

- 

0 

519.85 

591.68 

617.92 

640.62 

653.12 

671.95 






(9.73) 

(2.53) 

(1.63) 

(2.03) 

(2.19) 

(3.22) 

- 

- 

- 

- 

20 

473.01 

547.15 

572.38 

596.75 

609.85 

629.41 






(9.38) 

(2.54) 

(1.89) 

(1.85) 

(1.93) 

(2.65) 

- 

- 

- 

- 

40 

427.12 

500.85 

526.35 

549.89 

564.43 

585.13 






(7.87) 

(2.55) 

(1.83) 

(1.96) 

(2.26) 

(2.66) 

- 

- 

- 

- 

60 

378.13 

453.61 

480.29 

504.73 

518.80 

541.20 

553.63 

564.43 

572.46 

580.38 


(7.27) 

(2.89) 

(2.20) 

(2.27) 

(2.35) 

(2.67) 

(2.83) 

(3.32) 

(4.35) 

(5.31) 

80 





471.67 

495.48 

509.54 

519.92 

529.52 

537.30 


- 

- 

- 

- 

(2.62) 

(2.71) 

(3.24) 

(3.94) 

(3.93) 

(4.91) 

100 





425.38 

449.48 

466.34 

477.19 

486.61 

493.28 


- 

- 

- 

- 

(2.74) 

(2.92) 

(3.16) 

(3.18) 

(4.60) 

(5.36) 

120 





379.29 

403.94 

422.49 

432.59 

442.65 

450.24 


- 

- 

- 

- 

'(2.80) 

(2.95) 

(3.33) 

(3.51) 

(4.28) 

(5.06) 

140 





333.72 

357.35 

378.30 

389.32 

398.93 

406.40 


- 

- 

- 

- 

(2.78) 

(3.52) 

(2.97) 

(3.48) 

(4.29) 

(5.29) 

160 






312.08 

333.28 

344.88 

354.88 

362.48 


- 

- 

- 

- 

- 

(2.99) 

(2.92) 

(3.69) 

(4.52) 

(5.47) 

Standard deviations are given in parentheses. 
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Table 4. Average open-circuit voltages and standard deviations for 10 Sl-cm cells 


Tempera- 





Solar intensity,' 

1 mW/cm 2 





ture, °C 

5.00 

25.00 

50.00 

100.00 

140.00 

250.00 

400.00 

550.00 

700.00 

850.00 

-160 

748.38 











(95.76) 

- 

- 

- 

- 

- 

- 


- 

- 

-140 

775.15 

799.57 










(74.95) 

(21.75) 

- 

- 

- 

— 

- 

- 

- 

- 

-120 

761.75 

774.55 

771.32 









(50.28) 

(20.43) 

(21.03) 

- 

- 

— 

- 

- 

- 

- 

-100 

742.28 

752.01 

751.59 

748.30 








(30.28) 

(20.86) 

(21.84) 

(23.71) 

- 

- 

- 

- 

- 

- 

-80 

712.06 

729.50 

730.38 

727.95 

726.61 

703.65 






(18.27) 

(19.96) 

(22.26) 

(24.26) 

(23.24) 

(24.64) 

- 

- 

- 

- 

-60 

671.44 

701.81 

706.35 

707.47 

707.06 

682.18 






(12.73) 

(16.19) 

(20.01) 

(22.52) 

(23.47) 

(23.39) 

- 

- 

- 

- 

-40 

622.32 

664.17 

674.55 

681.40 

681.55 

652.53 






(10.86) 

(10.41) 

(15.09) 

(19.27) 

(20.88) 

(19.19) 

- 

- 

- 

- 

-20 

573.46 

620.71 

634.25 

646.38 

649.84 

617.85 






(9.47) 

(5.60) 

(9.30) 

(13.92) 

(15.97) 

(14.00) 

- 

- 

- 

- 

0 

522.92 

572.76 

590.65 

605.38 

610.96 

577.44 






(8.47) 

(3.12) 

(5.15) 

(8.61) 

(10.51) 

(9.29) 

- 

- 

- 

- 

20 

470.94 

524.32 

544.13 

560.22 

568.31 

533.95 






(7.45) 

(2.75) 

(3.39) 

(5.24) 

(6.45) 

(6.28) 

- 

- 

- 

- 

40 

419.75 

475.60 

495.21 

512.84 

522.12 

488.53 

497.22 

505.55 

510.02 

512.75 


(6.80) 

(2.99) 

(3.14) 

(3.72) 

(4.65) 

(4.97) 

(5.77) 

(6.87) 

(7.30) 

(7.47) 

60 

367.48 

425.50 

446.62 

466.42 

475.02 

440.91 

451.02 

459.19 

464.02 

467.15 


(5.97) 

(3.18) 

(3.45) 

(3.74) 

(3.95) 

(4.78) 

(5.25) 

(6.15) 

(6.39) 

(7.30) 

80 





427.74 

390.99 

403.42 

412.34 

417.72 

420.22 


— 

— 

- 

- 

(4.34) 

(4.88) 

(5.33) 

(6.16) 

(6.18) 

(6.99) 

100 





379.47 

341.22 

355.19 

364.48 

369.31 

372.91 


- 

- 

- 

- 

(4.38) 

(4.98) 

(5.38) 

(6.02) 

(6.66) 

(7.02) 

120 





330.45 

294.65 

308.23 

314.97 

321.32 

324.88 


- 

— 

- 

- 

(4.82) 

(5.03) 

(5.69) 

(6.03) 

(6.68) 

(7.03) 

140 





280.46 

244.52 

257.80 

266.90 

272.53 

277.72 


- 

- 

- 

- 

(5.11) 

(5.31) 

(5.91) 

(6.15) 

(6.55) 

(8.76) 

160 



— 

— 









_ 

— 

— 

“Standard deviations are given in parentheses. 
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The <ii and bi are presented in Table 5, The complete 
expression for V oc ( T,H ) is 

V 0C (T,H) = A(T) + B(T) log 10 H 

= a 0 + fli T + a 2 T 2 4- a 3 T 3 + o 4 T 4 + a 5 T 5 
+ (b 0 + b,T + b 2 T 2 + b 3 T 3 + bj* 

+ bsT 5 ) lo gl o H (10) 

Although the expression for V oc ( T , H) has many terms, 
it is well suited for computer computation. Thus, the 
entire array of parametric data for each cell type can be 
characterized by Eq. (10), together with the appropriate 
values for the 12 coefficients. Equation (10) was evalu- 
ated for both the 2 and 10 Q-cm cells at each of the 
temperature-intensity conditions under which experimen- 
tal data were taken. The percent difference between the 
computed V oc (using Eq. 10) and the actual V oc was cal- 
culated for each of the temperature-intensity conditions 
according to the expression 

percent difference (T *, H ,) 

_ calculated V oc (T u H ,) — actual V oc (T { , Hi) 

~ 100 1/2 calculated V 0C (T i ,H j ) + actual V 0C (T i ,H t ) 

( 11 ) 

The percent differences are tabulated for both 2 and 
10 Q-cm cells in Tables 6 and 7. 

It was found that, in general, the computed V oc agree 
quite well with the actual ones, the maximum difference 
occurring at 5 mW/cm 2 (1 to 6%). However, from 25 to 
850 mW/cm 2 , all of the percent differences are less than 
\% for 2 Q-cm cells, whereas all but four of the percent 
differences are less than 1% for the 10 Q-cm cells. Although 
Eq. (10) does not explicitly contain a temperature coeffi- 
cient term, a point value of fl can be determined at any 
temperature T 0 and intensity H 0 by differentiating Eq. (10) 
with respect to the temperature holding intensity constant 
and then evaluating the derivative at (T 0 , H 0 ). 

Also, Eq. (10) can be used to evaluate the average tem- 
perature coefficient /? (T, H) over a temperature range T, 
to T k at some intensity Hi. The appropriate expression is 

average temperature coefficient = /? (T, to T k , Hi) 


Again, this expression is well suited for computer 
computation. 

2. Development of l sc Model. The following expression 
was found to describe the I sc performance as a function of 
intensity at constant temperature 7\: 

I. c (Ti,H) = C T .H (13) 

Equation (13) fits the 2 Q-cm data quite well over the entire 
range of intensities (5 to 850 mW/cm 2 ), whereas the fit for 
the 10 Q-cm data was good only from 5 to 250 mW/cm 2 . 
However, since 10 Q-cm cells are not likely to be used for 
near-Sun missions, as discussed previously, they were not 
modeled for intensities greater than 250 mW/cm 2 . 

The temperature dependence of C (T) was found by the 
same method used for A (T) and B (T). The complete ex- 
pression for I ec takes the form 

I ec (T, H) = (c„ + c,T + c 2 T 2 + c 3 T 3 + c 4 T 4 + cj*) H 

(14) 

Although Eq. (14) does not explicitly contain a tempera- 
ture coefficient term, a point value of a can be determined 
at any temperature T 0 and intensity H 0 by differentiating 
Eq. (14) with respect to the temperature holding intensity 
constant and then evaluating the derivative at (T 0 ,H 0 ). 
Also, Eq. (14) can be used to determine the average short- 
circuit temperature coefficient a (T k — T,, H { ) over a tem- 
perature range T, to T lc at intensity Hi. The expression is 

s (T k - Tj, Hi) = [* (T k - T i) + c 2 (Tl- Tfj 
+ C3 (T 3 -T 3 ) + c 4 (r 4 -T 4 ) 

+ c 5 (T 3 -T 3 )]— (15) 

Equation (14) was used to calculate I 8C at each of the 
temperature-intensity conditions at which experimental 
data were taken, as shown in Tables 1 and 2. The percent 
differences between the calculated and actual I sc for both 
the 2 and 10 Q-cm cells were evaluated and are presented 
in Tables 8 and 9. It was found that most of the differ- 
ences fall in the range of 0 to 3% for both types of cells. 


_ ai (ft - Tj) + a 2 (Tj - TJ) + a 3 (Tj - T 3 ) + g 4 (T 4 - T*) + a „ (T 3 - T 3 ) 

T k - r, 

+ Ilk (T k ~ Tj) + b 2 (Tj ~ T)) + b 3 (T 3 - T 3 .) + b t (Tj - T 4 ) + b 5 (T 3 - Tp] log 10 H { 

T t ~ Tj 
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Table 5. Coefficients for analytical expressions b 


Voc expression 
V„c = A(T) + B(T)logioH 
A ( T ) — ao + aiT 1 + ( 12 T 2 + azT 3 + a 4 T 4 + asT 5 
B ( T ) = bo + biT 1 + b 2 T 2 + bsT 3 + b 4 T 4 + baT 3 


he expression 
he = C (T) H 

C (T) = co + ciT 1 + c 2 T 2 + C3T 3 + c 4 T 4 + c s T 3 


2 B-cm cells 0 


ao = 0.480653 + 03 

bo = 

0.801129 + 02 

CO = 

0.914727 + 00 

ai = -0.241702 + 01 

bi = 

0.118309 + 00 

Cl = 

0.108713 - 02 

a 2 = 0.294965 - 02 

&2 = 

-0.224011 - 02 

C2 = 

-0.695706 - 05 

as = -0.206159 - 04 

bz — 

0.140478 - 04 

‘ C3 = 

0.226603 - 07 

tf4 = -0.184560 - 06 

b 4 = 

0.816163 - 07 

C4 = 

0.171090 - 09 

as = 0.114636 - 08 

bs = 

-0.548025 - 09 

C5 = 

-0.144039 - 11 

10 Jl-cm cells' 1 

do = 0.507642 + 03 

bo = 

0.480243 + 02 

CO = 

0.100584 + 01 

ai = -0.300212 + 01 

bi = 

0.455891 + 00 

Cl = 

0.143787 — 02 

a 2 = 0.407518 - 02 

bi — 

-0.470785 - 02 

C2 = 

-0.315471 - 05 

as = 0.459916 - 04 

bz = 

-0.909273 - 05 

C3 = 

-0.121720 - 06 

a 4 = -0.323551 - 06 

b 4 = 

0.205860 - 06 

C4 = 

0.538927 - 10 

as = -0.164756 - 09 

bs = 

-0.323078 - 09 

C5 = 

0.462234 - 11 


a Coefficients are in “E” format. The two digits on the right, including the sign, are powers of 10 which multiply the six digits to the left 
of the powers. 

b Units: V oc in mV, he in mA, H in mW/cm 2 , T in °C. 

c See Table 6 for applicable temperature intensity range for V oc , Table 8 for he- 
d See Table 7 for applicable temperature intensity range for V 0 c, Table 9 for he. 


By comparing Tables 6 and 7 with Tables 8 and 9, it can 
be seen that the analytical expressions fit the V oc data 
better than the I sc data. This is to be expected, since tem- 
perature is more accurately controlled than intensity dur- 
ing the parametric cell testing; I sc depends linearly on 
intensity, while V oc is a logarithmic function of intensity. 
For example, at 140 mW/cm 2 , a ±3% error in intensity 
causes a corresponding error in l sc , while the effect would 
be much less than 1 % on V oc . 

Method 2 was used together with prelaunch trajectory 
and temperature predictions to predict the Mariner 9 
I BC — V oc transducer performance as a function of helio- 
centric distance. Data from the I sc — V oc transducer dur- 
ing spacecraft flight, together with accurate transducer 
predictions, have proven to be of significant value in pro- 
viding engineering information on solar array performance 
characteristics. 


C. Method 3: Computer Program Interpolation Scheme 

This method was developed in order to provide a tech- 
nique for predicting the solar electrical performance over 
a wide range of temperature-intensity conditions and to 

(1) Preserve the original identity of the parametric data 
by minimizing the use of derived quantities. 


(2) Utilize only those parametric data which cover the 
range of the predictions being made. For example, 
for a Mercury mission, only the parametric data in 
the intensity range of 140 to 850 mW/cm 2 would 
be used. 

(3) Provide a method for propagating parametric data 
uncertainty (as reflected by standard deviations) 
through prediction calculations and give some esti- 
mate of the uncertainty of the prediction. 

This method uses a computer program interpolation 
scheme. The A T . and B T . are found for the V oc data and 
C Ti are found for the I sc data in exactly the same fashion 
used in method 2. However, only the parametric data 
which cover the range of the predictions being made are 
utilized. Once the A Ti , B T{ , and C Ti are determined at each 
parametric temperature, Eqs. (6) and (13) are evaluated 
for each set of A, B, and C (one set per parametric tem- 
perature) at the prediction intensity. This method differs 
from method 2 in that no attempt is made to find the 
coefficient’s functional dependence on temperature. The 
temperature dependence of V oc and I ec at the prediction 
intensity is found directly by a “French Curve Fit” (seg- 
mented cubic polynomials) of the results obtained by 
evaluating Eqs. (6) and (13), respectively, at each of the 
parametric temperatures. The French Curve is then evalu- 
ated at the prediction temperature. 
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Table 6. Comparison of computed and actual V oc , 2 si-cm cells (% difference) 


Temperature, 

°C 




Solar intensity, mW/cm 2 





5.00 

25.00 

50.00 

100.00 

140.00 

250.00 

400.00 

550.00 

700.00 

850.00 

-160 

-1.00 



— 

— 

— 

— 

— 

— 

— 

— 

-140 

3.57 

-0.83 

— 

— 

— 

— 

— 

— 

- 

- 

-120 

4.76 

0.02 

0.01 

— 

— 

— 

— 

— 

- 

- 

-100 

5.20 

0.03 

-0.55 

0.33 

— 

— 

— 

— 

- 

- 

-80 

5.07 

0.38 

-0.15 

0.06 

0.02 

— 

— 

— 

— 

- 

-60 

4.28 

0.12 

-0.18 

-0.13 

0.22 

0.36 

— 

— 

- 

- 

-40 

3.80 

0.09 

-0.30 

0.11 

-0.23 

0.17 

— 

— 

— 

- 

-20 

3.19 

0.03 

-0.19 

0.24 

-0.16 

0.10 

— 

— 

— 

- 

0 

3.18 

0.16 

-0.19 

0.04 

-0.08 

0.12 

— 

— 

— 

- 

20 

3.61 

0.07 

-0.05 

-0.01 

-0.20 

-0.03 

— 

— 

- 

- 

40 

4.00 

0.25 

0.10 

0.32 

-0.13 

-0.12 

— 

— 

— 

— 

60 

5.17 

0.54 

0.15 

0.24 

-0.15 

-0.43 

0.39 

0.50 

0.60 

0.43 

80 

— 



_ 

— ' 

0.08 

-0.46 

0.16 

0.39 

0.23 

0.10 

100 







— 

0.14 

-0.39 

-0.22 

0.01 

-0.07 

0.06 

120 

_ 







0.20 

-0.33 

-0.38 

0.16 

0.01 

-0.00 

140 







_ 

0.17 

0.11 

-0.41 

0.09 

0.12 

0.22 

160 

— 





— 

— 

0.19 

-0.35 

0.12 

0.11 

0.23 



Table 7. Comparison of computed and actual V 0 c, 

10 Sl-cm cells (% difference) 



Temperature, 

°C 




Solar intensity, mW/cm 2 





5.00 

25.00 

50.00 

100.00 

140.00 

250.00 

400.00 

550.00 

700.00 

850.00 

-160 

0.42 





— 

_ 

- 

— 

— 

— 

— 

-140 

0.04 

-1.86 

— 

- 

— 

— 

— 

— 

— 

- 

-120 

2.17 

-0.22 

-0.11 


— 

— 

— 

— 

— 

— 

-100 

2.85 

0.46 

0.05 

0.02 

— 

— 

— 

— 

— 

— 

-80 

3.11 

0.27 

-0.04 

0.11 

0.21 


— 

— 

— 

— 

-60 

3.37 

-0.19 

-0.47 

-0.26 

-0.02 

0.77 

— 

— 

— 

— 

-40 

3.86 

-0.15 

-0.64 

-0.61 

-0.12 

0.65 

— 

— 

— 

— 

-20 

3.66 

0.01 

-0.36 

-0.50 

-0.19 

0.82 

— 

— 

— 

— 

0 

3.44 

0.35 

-0.24 

-0.28 

-0.04 

0.80 

— 

— 

— 

— 

20 

3.59 

0.46 

-0.14 

-0.04 

-0.05 

0.78 

— 

— 

— 

— 

40 

4.06 

0.53 

0.12 

0.12 

-0.02 

0.53 

— 

— 

— 

— 

60 

5.68 

1.08 

0.27 

-0.18 

-0.18 

0.09 

0.76 

0.71 

1.03 

1.46 

80 

— 

— 

— 

— 

-0.45 

-0.16 

0.18 

0.12 

0.36 

0.72 

100 

— 

— 

— 

— 

-0.45 

0.15 

-0.16 

-0.49 

-0.39 

0.11 

120 

— 

— 

_ 

— 

-0.14 

0.68 

-0.18 

-0.69 

-0.46 

-0.20 

140 

— 

— 

_ 

_ 

. 0.25 

0.24 

-0.45 

-0.11 

-0.25 

0.12 

160 

— 

_ 

_ 



— 

0.09 

0.02 

-0.07 

0.34 

0.41 
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Table 8. Comparison of computed and actual l„, 2 n-cm cells (% difference) 


Temperature, 

°C 




Solar intensity, mW/cm 2 





5.00 

25.00 

50.00 

100.00 

140.00 

250.00 

400.00 

550.00 

700.00 

850.00 

-160 

1.24 

— 

— 

— 

— 

— 

— 

— 

— 

— 

-140 

-3.72 

-0.23 

— 

— 

— 

— 

— 

— 

— 

— 

-120 

1.44 

0.20 

-2.42 

— 

— 

— 

— 

— 

— 

— 

-100 

2.90 

1.85 

-0.98 

-1.00 

— 

— 

— 

— 

— 

— 

-80 

3.82 

2.08 

1.07 

1.49 

0.05 

— 

— 

— 

— 

— 

- -60 

4.08 

3.78 

1.63 

0.44 

-0.79 

-1.70 

— 

_ 

— 

— 

-40 

3.95 

3.27 

1.89 

0.48 

-1.39 

-2.89 

— 

— 

— 

— 

-20 

2.74 

4.60 

3.06 

0.80 

-1.15 

-2.63 

— 

— 

— 

_ 

0 

3.64 

2.11 

4.07 

-0.53 

-0.99 

-2.45 

— 

— 

— 


20 

4.59 

0.50 

2.83 

-0.07 

-1.78 

-2.42 

— 

— 

— 

— 

40 

-2.41 

0.72 

1.98 

0.06 

-1.51 

-0.80 

— 

— 

— 

_ 

60 

-1.97 

0.80 

1.84 

0.33 

0.27 

-0.80 

2.62 

2.73 

0.05 

-2.25 

80 

— 

— 

— 

— 

-0.25 

-1.28 

2.78 

2.34 

1.03 

-2.10 

100 

— 

— 

— 

— 

0.79 

-1.73 

3.44 

2.43 

1.37 

-2.10 

120 

— 

— 

— 

— 

0.33 

-1.48 

1.63 

2.68 

0.76 

-2.35 

140 

— 

— 

— 

— 

-1.29 

-2.94 

-0.69 

0.33 

1.09 

-2.75 

160 

— 

— 

■ — 

— 

— 

-5.68 

-3.16 

-0.41 

-1.54 

-4.80 



Table 9. Comparison of computed and actual /«„, 

10 O-cm cells (% difference) 



Temperature, 

°C 




Solar intensity, mW/cm 2 





5.00 

25.00 

50.00 

100.00 

140.00 

250.00 

400.00 

550.00 

700.00 

850.00 

-160 

-16.27 

— 

— 

— 

— 

— 

— 

— 

_ 



-140 

-5.41 

-1.26 

— 

— 

— 

— 

— 

— 

— 



-120 

-3.26 

1.44 

0.94 

— 

— 

— 

— 

— 





-100 

-3.73 

1.17 

0.59 

0.10 

— 

— 

— 

— 

_ 



-80 

-4.05 

0.52 

0.02 

-1.09 

-1.43 

— 

— 

— 

— 



-60 

-4.45 

-0.98 

0.63 

-2.80 

-1.69 

-1.91 

— 

— 

— 



-40 

0.50 

-0.55 

0.44 

-1.92 

-1.70 

-1.34 

— 

— 

— 



-20 

3.20 

0.25 

0.70 

-1.17 

-0.48 

-0.55 

— 

— 

— 



0 

4.46 

1.43 

3.03 

0.03 

-0.16 

-0.61 

— 

— 

— 



20 

7.47 

2.15 

2.33 

0.58 

-0.71 

-0.17 

— 

— 

— 

— 

40 

7.21 

2.46 

2.11 

1.02 

1.21 

0.34 

— 

— 

— 

— 

60 

6.11 

2.10 

1.31 

0.54 

0.82 

-0.12 

2.58 

1.20 

3.21 

5.69 

80 

- 

— 

— 

— 

-0.36 

-1.50 

1.62 

0.53 

3.02 

6.53 

100 

— 

— 

— 

— 

-2.11 

-1.16 

0.72 

0.05 

3.53 

8.21 

120 

— 

— 

— 

— 

-4.91 

-1.77 

0.95 

1.60 

6.75 

13.34 

140 

— 

— 

— 

— 

-2.39 

-1.77 

5.02 

9.90 

16.33 

23.78 

160 

- 

- 

- 

- 

- 

13.26 

19.39 

26.32 

34.67 

43.01 
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Parametric data uncertainty is propagated through the 
prediction calculation by means of a series of repeated 
predictions, as described in the preceding paragraph. 
However, between predictions, each parametric data point 
is modified by adding to it the product of its standard 
deviation and a restricted random number function R, 
with an average value of zero and a standard deviation 
of one. After an appropriate number of iterations, the 
average and standard deviation of the various predictions 
are taken. 

This method is the most accurate of the three; however, 
it is also the most inconvenient to use as the entire array 
of parametric data characterizes a particular cell type. 

V. Conclusions 

The following conclusions are based on the results of 
the experimental solar test program and the solar cell pre- 
diction models discussed in this report: 

(1) The 2 Q-cm cells are suitable for standard cells up 
to intensities of at least 850 mW/cm 2 (at a reference 


temperature of 28°C); 10 Q-cm cells are suitable up 
to 250 mW/cm 2 . Beyond that, the I ec of the 10 Q-cm 
cells falls off from a linear function of intensity. 

(2) For near-Sun missions, 2 Q-cm cells should be 
selected for use over 10 Q-cm cells, since above 
250 mW/cm 2 the P mp of the latter falls off with 
respect to the P mp of the former. 

(3) Deviation in P mp at the high- and low-intensity 
extremes (as reflected by standard deviations) indi- 
cates that for inner and outer planet missions fur- 
ther attention must be focused on cell fabrication 
processes, proper selection of materials, and cell 
screening techniques in order to ensure optimum 
solar array performance. 

(4) Analytical expressions derived from parametric data 
(method 2) and a computer program interpolation 
scheme (method 3) can be used to accurately repro- 
duce the laboratory solar cell parametric data over 
a wide range of temperature and intensity. 
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